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Abstract 



The probability distribution dw/df as a function of the ratio / = no/ntot of 
\ the neutral to total multiplicities is calculated for the classical pion fields quickly 

varying in space and time. 

CO 

1. It is widely discussed in the recent years that the production of a classical pion 
\ field leads to the ratio 

>: / = -=*■ (i) 

Q-T 

of the multiplicity of the neutral pions to the total multiplicity which completely differs 
from what one expects for the usual mechanisms of the pion production. For the large 
values of n tot the standard mechanisms predicts the probability 

while for a particular classical state known as the "disoriented chiral condensate" (DCC), 
the distribution is 

— - — (3) 

(see, for example, Q and a mini- review || and the references therein). 

In this letter we shall discuss a class of classical chiral fields described in and show 
that the distribution (3) is correct only for the fields slowly varying in space and time (for 
instance, for the DCC which is an exact constant) whereas for the fields which vary rapidly 
enough the predicted distribution is quite different (the plot of the latter distribution is 
shown in Fig.l). 

The classical solutions of the non-linear a-model found in [[3] can be represented in 
the following form. As usually the four fields a, tt are constrained by the relation: 

- 2 + - 2 = /. 2 - (4) 
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If we denote $ a = f n 1 (n,a), a = 1, 2, 3, 4, the solution is 

$ a = A a cos(6(x) +<f a ) , (5) 
where the amplitudes A a and the phases <p a obey the equations: 

E^a = 2, £ ^e 2 ^ = , (6) 

a=l a=l 

while the function #(x) satisfies the free wave equation: 

d 2 6{x) = . (7) 

The derivation and the discussion of these solutions see in ref.0. 
Let us introduce the complex-4- vector in the chiral 0(4) space: 

Z a = A J"* = X a + tY a . (8) 

Eqs. (6) can be rewritten 

X 2 a = l, Y a 2 = l, X a Y a = 0. (9) 

(HereX a 2 = ELi^, etc.). 

The 0(4) invariant dynamics for the production of the states characterized by the pa- 
rameters X a , Y a can lead to the distribution function which depends only on the invariants 
X li Y a ,X a ■ Y a : 

dw = p{Xl Yl X a ■ Y a ) d A Xd 4 Y . (10) 

Due to Eqs. (9) the function p is proportional to the 5-functions: S(X 2 — 1)5(Y 2 — 
l)5(X a Y a ), and since all the invariants turn out to be fixed (10) takes the form: 

dw = const 5(Xl - 1)5(Y 2 - l)5(X a Y a )d A Xd A Y . (11) 

The number of the neutral pions no and the total number of the pions n to t are the 
functions of X a , Y a . The distribution in the variable /, the fraction of neutral pions, is 
given by: 

^ = const fs(f- 2^ X Y Y) ) 5{Xl ~ ~ mXaYa)*X*Y. (12) 

The ratio n /n to t is proportional to ratio of the squared amplitudes of the correspond- 
ing pion fields integrated over the production volume: 

n (X,Y) _ Jd 3 r$ 2 



n to t(X,Y) /rf 3 r[$ 2 + $i + $i 
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According to (5): 

$ 2 = A 2 cos 2 (d{x) + ^) , 2 = 1,2,3. (14) 
Two different limiting cases are possible. 

If 6{x) varies slowly through the production volume, 6{x) ~const, one can redefine the 
phases and put 9 = 0. Then 

$ 2 = A 2 cos 2 ^ = X 2 , 



n (X,Y) X. 



(15) 



n tot (X,Y) " XI + XI + XI ■ 

This case of the constant field is actually the case of the DCC. 

Substituting (15) into Eq.(12) we can easily perform an integration over d 4 Y and the 
three of the four d 4 X integrations to get: 

dw f +1 1 

— = const / d(cos9)5(f — (cos#) 2 ) = const —=. (16) 

df J-i V/ 

Here 

cos 9 = , 2 X \ = • (17) 
y X 1 + X 2 + X 3 

Thus we see that the inverse square root low comes out when the field varies slowly 
through the production volume. 

For the opposite case we write cos 2 (9(x) = 1/2(1 + cos(26>(:r) +2</?j)) and assume 
that the term linear in cos(26 l + 2</?j) is integrated out in Eq.(13). Then 

n (X, Y) _ A 2 _ XI + Y 2 ^ 2 = j, x ^ y 2 = Y, ■ (18) 



n to t(X,Y) A 2 X 2 + Y 2 



i=l i=l 



The integral, which appears when this expression is substituted into Eq.(12), is less 
trivial then the one for the previous case. In fact it appears possible to perform the 7 of 
the 8 integrations d 4 Xd 4 Y analytically. This is described in the second part of this letter. 
A reader who is not interested in the technical details turn to the plot of J(f) = dw/df 
shown in Fig.l. Note also that 



7T 

tan- 



+ \ = 2.058 , J(l) = . (19) 



3 ' \2J 2y/2 

The discontinuity of dJ(f)/df at / = 1/2 is elucidated qualitatively in the next section 
of the letter. 

2. The properly normalized distribution which emerges from Eqs. (12) and (8) is 
f = J(f) = ^Jd 4 Xd 4 Y5^f-^^y(x 2 a -l)5(Y a 2 -l)5(X a Y a ), 

f Af)df = 1. (20) 
Jo 



We split d A X = d 2 x±d 2 x\\ and d A Y = d 2 y±d 2 y\\, where x± = (Xi,X 2 ), x\\ = (X^X^), 
y± = (Yi,Y 2 ), y\\ = (13,14), and perform the integrations cPx±cPy± using the three 5- 
functions. One easily obtains: 



J(f) = ^.{^(f-J^y?) 



X 



- - y\) - (%1) 2 ]Q(1 - xg)8(l - yjf) 

where O's are the step-functions. The last two O functions can be rewritten in the 
following way. To have 1 — x 2 > and 1 — y 2 > it is necessary and sufficient to have 
the product of these factors and their sum to be positive. The product (1 — x 2 )(l — yjj) is 
anyway positive due to the first © function. Therefore we can change 0(1— x|)0(l — y 2 ) — > 
9(2 yjf). 

Using Eq.(8) we rewrite the integrations d 2 x\\d 2 y\\ = AsdAsdipzA^dA^dip^ and readily 
get for J(/): 

J U) = i r"^ r dX(l - A) e[l-2A-2/(l-A) + 4A(l-A)/ 3m Vl^ 
nJo Jo yjl - 2A - 2/(1 - A) + 4/A(l - A) am 2 <p 

In Eq.(22) we used \ = 1/2 A 2 and tp — tp^ — tp^. 

Two limiting cases are immediately obtained from (22) . For / = 1 the argument 

of the O function is negative (except for one point: A = 1/2, sin 2 ip = 1) and therefore 
J(l) = 0. For / = one has: 

V ; Jo VT^2X 3 V ; 



For the arbitrary / it is convenient to perform the A integration to get: 

n /2 A=A max , or 1 

T/ , N 4 r , f VaA 2 + b\ + c 1 + b/2a 2a\ + 6 

J(j) = — / a</? < arcsm 



A=A min , or 



Here 



71 J [ -a ^f^a ^Jb 2 - 4ac 

-Af sin 2 tp, b = -2(1 - /) + 4/ sin 2 ip, c = 1 - 2f , 



(24) 



b b 2 -4ac 
A max , min - -- ± \/^^- (25) 

The limits (A m i n , A max ) in Eq.(24) should be understood in the following way. If b 2 — 
4ac < then J = (since the argument of the 0(aA 2 + 6A + c) in Eq.(22) is negative). For 
b 2 — 4ac > the integration in A is taken between A m j n and A max only if these quantities 
are inside the interval (0,1). 
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The last integration in (24) over ip has been performed numerically and led to the plot 
for J(f) depicted in Fig.l. We would like to make two comments concerning this plot. 

First, one can analytically calculate </(/) for / = 1/2. From Eq.(24) one obtains: 



,7(1/2) = V2 J d<p 



w/2 

1 1 

+ 



tt/4 



sin (p 2 sin ip 



In tan - + - = 2.058 . (26) 



2y/2 8 2 



Second, to understand qualitatively the discontinuity of dJ(f)/df at / = 1/2 one can 
notice that for the fixed value ip = 7r/4 the integrand in Eq.(24) has a gap for / = 1/2 + e 
when e — > +0. 

Indeed, for ip = 7r/4 one has a = — 1 — 2e, b = +4e, c = —2e. Since it will be 
immediately evident that only small A ~ are essential the polynom aX 2 + bX + c ~ 
—A 2 — 2e. The integral (22) in A (i.e. the integrand in(24)) is then 

l/'dA 9 '-^-^ = I**? £ <° (27) 
7T Jo V-2e - A 2 \ for e > . K ' 

The integration over (p (within the interval Acp ~ y/i, where bX < \c\ smears out the 
discontinuity in J(f) but the derivative dJ/df still has "singularity" dJ/df ~ —Q(e)/\/e 
at e -»• +0. 

3. In ref. we discussed the possible production mechanisms and the signatures which 
can be used for the observation of the classical pion field. One of the conclusions of this 
paper was that it is not easy to observe the DCC, or almost constant pion field, due to a 
simple reason: the energy density for such a state is so small that no more than tenths of 
the pions can be produced through the DCC decay. We claimed that the observation of 
the field varying in space and time may be more promissing since in this case the energy 
density can be much higher. We see now that the distribution in / = n /n tot for such 
fields may be quite different from what is usually supposed, i.e. from dw/df ~ l/\/J, 
and that for a class of the fields described in dw/df is given by Eq.(24) (Fig.l). It also 
seems that the latter fields are closer to those discussed in [[|] that the constant DCC.Q 
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1 We are grateful to K.Rajagopal for this comment. 
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Figure 1: The probability distribution dw/df over the neutral to total multiplicities ratio 
/ = n /n tot for classical pion fields quickly varying in space and time. 
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